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ABSTRACT 

Previous studies suggest that the planet-forming disks around very-low-mass stars/brown 



dwarfs may be flatter than those around more massive stars, in contrast to model predictions of 
t — ■ ■ larger scale heights for gas-disks around lower-mass stars. We conducted a statistically robust 

study to determine whether there is evidence for stellar-mass-dependent disk structure in pianet- 
te — ■ forming disks. We find a statistically significant difference in the Spitzer/IRAC color distributions 
of disks around very-low-mass and low-mass stars all belonging to the same star-forming region, 
the Chamaeleon I star-forming region. We show that self consistently calculated disk models 
cannot fit the median spectral energy distributions (SEDs) of the two groups. These SEDs can 
be only explained by flatter disk models, consistent with the effect of dust settling in disks. We 
find that relative to the disk structure predicted for flared disks the required reduction in disk 
scale height is anti-correlated with the stellar mass, i.e. disks around lower-mass stars are flatter. 
Our results show that the initial and boundary conditions of planet formation are stellar-mass- 
dependent, an important finding that must be considered in planet formation models. 

Subject headings: protoplanetary disks — planets and satellites: formation — infrared: stars — stars: 
low-mass — brown dwarfs 

1. Introduction grains into planetesimals and for the assembly of 

planets from planetesi mals (e.g. see reviews in 



Circumstellar disks around young stars provide lApai fc Lauretta] l2010l) . Cool stars, with spectral 



the raw material for planet formation. In addi- types M0 or laterj account for more than 80% of 

tion, the evolution of protoplanetary disks sets the the galactic stellar population. Their abundance 

boundary conditions for the agglomeration of dust _ comb i ne d with the suite o f recent disco veries 

of planets around M-dwarfs ( Gaudi et al.l 12008) 
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hosts in the Galaxy. But how do these planetary 
systems differ from those around sun-like stars? 

Lacking observational data, some models of 
planet formation assumed that disk properties 
are independent o f the stellar properties (e.g. 



Kornet et al.l 120061 ) . Disk mass measurements, 



however, revealed that the disk/stellar mass ratio 

is si milar across a broad range of stellar masses 

e.g.lKlein et al.ll2003l; lAndrews fe Williams! I200J 



Scholz et al.l 120061: IBouv et al.l l2008J h This led 
to refined models that explained the observed 
low occur rence rate of gas gia n t planets around 
M-stars (lLaughlin et all 120041: llda fe Linl 120051: 



Johnson et al.l2007tlBenz et all2008l : lKennedv fe Kenvon 
20081 ). However, all other disk parameters are 
currently assumed to be essentially stellar-mass- 
independent. 

Recent observational studies, however, revealed 
stcllar-mass-dependence for most disk proper- 
ties. The typical lifetime of optically thick 
dust disks around cool stars and brown dwarfs 
seems to be at least a factor of two or three 
longer dSterzik et al.lliool ICarpenter et alJl2006t 
Riaz et al.l 20061) than aroun d sun-like stars (see 
Pascucci fe Tachibanall2010l for a review). Ev- 
idence mounts that stellar mass accretion rates 
scale as M ex M 2 , resulting in very-low-mass ac- 
cretion rates at the low-end of the substellar mass 
spectrum (e.g. Natta et al. l2004t Mohantv et al 



20051: iMuzerolle et al.l .2005: H erczeg et aTT l2009) 



Surprisingly, dust processing in the mid-infrared- 
emitting regions of brown dwarfs and low-mass 
stars is more advanced than in disks around sun- 

20021 



like stars of similar 



age 



(e.g. Apai et al 



2004, 2005; Kcssler-Silacci et al. 2006; Bouv et al 
2008; IPascucci et al.ll2009t iRiaz et a l. 2009). Ex- 



citingly, Spitzer gas-line spectroscopy revealed a 
prominent difference in the organic chemistry in 
disks around very co ol stars and sun-like stars 



(jPascucci et al.1 l2009h . In addition, several ear- 



lier studies have suggested that disks around cool 
stars and brown dwarfs may be flatter than those 
around more massiv e stars (e.g. lApai et al.l 12002 ; 
Pascucci et al.l l2003t lApai et al.l 120051) . although 



disk models predict the opposite (jWalker et al 
20041 ). 



The goal of this work is to provide a statistically 
robust evaluation of the hypothesis that the low- 
est mass stars have, on average, flatter disks than 
their higher mass counterparts. First, we present 



results from a statistical study of Spitzer/IRAC 
colors of the complete census of the Chamaeleon I 
star-forming region. This comparison reveals a 
prominent difference between the color distribu- 
tion of the very-low- and low-mass stars in the 
sample. Next, we use a grid of simple disk models 
to demonstrate that the difference in the distribu- 
tions of IRAC colors is well explained by the lower 
occurrence rate of flared disks around the lowest- 
mass stars. These findings further reinforce the 
stellar-mass-dependent evolution of protoplane- 
tary disks, a key — and yet unaccounted for - 
factor to consider in planet formation models. 

2. Data Analysis 

In this study we focused on the Chamaeleon I 
star-forming re gion (Cha I, ^200 known members 
Luhmanl l2008al lbl) . which provides an ideal set of 
targets. The membership of the cluster can be 
considered virtually complete for masses greater 
than ~16 Mj, providing an unbiased sample of 
roughly coeval, well-characterized stars and disks 
( Luhman 20071: Apai et al. 2005t Pascucci et al 



20081 120091 iNatta et all 120041: IRiaz et al.l 12009) 
While the mean age of the very-low- and low- 
mass stars may be slightly different in Cha I, the 
difference is thought to be lowe r than the esti- 
mated age spreads for each group ( Luhmanll2008cf ) 
and therefore can be assumed coeval for the pur- 
poses of this study. 

We used phot ometric data from the Spitzer 
Space Telescope ((Wernerl 120051 ) obtained by the 
Infrared Array Camera (IRAC. Fazio et al.ll2004f ) 
at 3.6, 4.5, 5.8, and 8 /im wavelengths and 
the Mu ltiband Ima ging Photometer for Spitzer 
(MIPS, iRiekd 12004 at 24 /xm. The IRAC and 
MIPS observations were performed within a ra- 
dius of 3° from a = ll' l 07 m 00 s , S = -77°10'00". 
We complemented these data with JHK S mea- 
surement s from the Two Micr on All Sky Survey 
(2MASS, ISkrutskie et~ai1l2006l). W e use the pho- 
tometric data fromlLuhmanl (|2008alfbh and 2MASS 
(|Skrutskie et al.ll2006l) with the p hysical parame- 
ters of the members determined bv lLuhmanl ( 2007 . 
2008bJ). For each source we adopted spec troscopi- 
cally de termined J-band extinction from iLuhmanl 
(|2008al lbl) and dereddened the measured magni- 
tudes using the extinction law derived by the 
linear interpolation of A(X)/A(J) from iMathis 
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Fig. 1. — The distribution of the IRAC-IRAC colors. The solid lines arc the distributions of the low-mass 
stars, while the dashed lines represent the VLMOs. The peak close to mag represents the stars without 
infrared excess. The distributions of the disk-bearing VLMOs are consistently shifted towards bluer colors 
on each diagram. 



( 19901 ). For the purpose of comparing spectral 
energy distributions we calculated the flux densi- 
ties in each band using zero points o f the IRAC 
and MIPS instruments ( Cohen et al.l 1200 3: Ricke 
20041 : iReachet al.l l2005h . The spectral types of 
the stars and brown dwarfs in our sample range 
from B6.5 to M9.75 (|Luhmanll2007h . We divided 
these objects into two groups: A) low-mass stars 
(from G5 to M4.5, 110 stars); and, B) very-low- 
mass objects (VLMOs, from M4.75 to M9.5, 90 
stars). The separation between the two groups 
corresponds to T e g = 3, 170^. 

3. Results 

We computed the IRAC colors for all objects 
as the difference between a magnitude at a short 
and at a longer wavelength. Fig. [1] shows the 
IRAC-IRAC color distributions of the low-mass 
stars and VLMOs. Objects with colors close to 
zero are consistent with pure photospheric emis- 
sion, while objects with red colors (higher flux at 
longer wavelengths) have excess emi s sion a rising 
from a dust disk. lHartmann et al.l (|2005>) have 
evaluated IRAC-IRAC colors as disk indicators 
and have shown that the [5.8] — [8] > 0.35 cri- 
terion is efficient in discriminating between stars 
with and without disks. We adopt this criterion 
for our sample of Cha I objects leading to 78 disk- 
bearing and 90 diskless stars. There are 32 stars 
that do not have either 5.8 or 8 /mi measurements, 
and thus cannot be classified. 



We find that disks in the VLMO group have 
bluer IRAC-IRAC colors than disks in the low- 
mass sample (see Fig. [I]). To verify the statisti- 
cal significance of this diff erence we performed a 
Kolmogorov-Smirnov test (jVetterling fc Flannerv 
2002h on the 47 low— mass and 31 VLMO disk- 



bearing members of the cluster. 

In the case of the IRAC-IRAC color distribu- 
tions the probability that the two samples are 
drawn from the same parent population are be- 
low 0.01 (see Tabled]). The low probabilities in 
all colors provide a solid proof for a different color 
distributions in the two groups. Furthermore, we 
also performed a K-S test for the IRAC-MIPS 
color distributions. For the IRAC-MIPS distribu- 
tions the probability of drawing the two distribu- 
tions from the same parent population is between 
0.1114 and 0.0174, suggesting a possible difference, 
but not reaching the satisfactory confidence level. 

When interpreting this result, however, an im- 
portant fact has to be kept in mind. While the 
IRAC maps are sensitive to the photospheres of all 
stars and brown dwarfs in Cha I, the MIPS mea- 
surements are sensitivity-limited and are therefore 
biased. The detection limit of the MIPS measure- 
ments introduces a bias for the fainter VLMOs: 
the fainter, less flared disks and/or diskless ob- 
jects are often missing from the MIPS maps. The 
absence of the stars with less flared disks skews 
the IRAC-MIPS color distribution toward flared 
disks and reduces the observable difference in the 



Table 1 
Probabilities that the IRAC colors of the disk bearing VLMOs and low-mass stars are 

DRAWN FROM THE SAME PARENT POPULATION, AS DETERMINED WITH A KOLMOGOROV-SmIRNOV TEST. 



IRAC Colors [3.6]- 



[4.5]-[8] [5.8]-[8] 



Probability 0.00219 0.00023 0.00074 



IRAC-MIPS distributions of the samples. In spite 
of the biased MIPS photometry the K-S test still 
suggests a difference between the two samples. 

We also find a difference in the IRAC-IRAC 
color distributions of the diskless stars. The color 
distributions of the diskless VLMOs are shifted 
towards redder colors, most likely due to strong 
molecular absorption bands in the atmospheres of 
the very-low-mass stars. Note, however, that this 
shift is to the opposite direction to the one seen in 
the disk-bearing sample. 

There are three stars in our sample with 
IRAC colors consistent with pure photospheres 
(-0.0 mag), but with large [8]- [24] color (>2 
mag). These so urces could be sta rs with tran- 
sition disks (e.g. IStrom et al.l Il989l for the def- 
inition). These stars do not alter our statistics 
significantly, but they are interesting for further 
studies. The 2MASS identifiers of these stars 
are J11022491-7733357, J11071330-7743498, and 
J11124268-7722230. 

In summary, the K-S test convincingly demon- 
strates that the color distribution of the disks 
around the VLMOs and low-mass stars is dif- 
ferent: VLMOs have disks with distinctly bluer 
IRAC colors. Earlier studies of small disk samples 
have hinted on such differences, often thought to 
be linked to the effects of dust settling on disks 
( Apai et al.ll2004L[2 005: Ke ssler-Silacci et a 1120061 



Dullemond fc Dominikll2004ri iFurlan et alJI^Ooi T 
Seeking an explanation for the observed color dif- 
ference, in the next section we will use a grid of 
simple disk models to verify this explanation and 
explore whether other possibilities are also viable. 

4. Disk Model Grid 

Fig. [5] shows the spectral energy distributions 
(SEDs) of all disk-bearing members of the low- 
mass and VLMO groups, normalized to their H- 
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Fig. 2. — Spectral Energy Distributions of the 
disk bearing members of the VLMO and low-mass 
samples normalized to the H band. 



band flux density. The 24 (im flux densities are not 
shown on this figure, because of the biased subset 
of the VLMO groups, as discussed in Section 3. 

While the median SEDs of the two groups have 
similar shape, the mean SEDs are different. In 
the case of the VLMOs the medians and means 
of the fluxes have similar values at a given wave- 
length, while in the case of the low-mass sample 
the mean IRAC values are higher than the me- 
dian IRAC values. This indicates more significant 
broadening of the flux density distributions to- 
wards higher fluxes than in the case of the VLMOs 
(see Fig. [2]). The higher broadening of flux density 
distributions of the low-mass sample is consistent 
with the results of the K-S tests: the flared disk 
to flatter disk ratio is higher in the case of the 
low-mass stars than for VLMOs. 



For modeling the disks we will use the median 
SED of the unbiased low-mass sample (including 
the 24 /im wavelength) and the median SED of the 
unbiased VLMO sample (JHK and IRAC bands), 
with the median of the biased 24 ^m flux density 
as an upper limit. The median 24 /jm flux den- 
sity of the sensitivity-limited (and therefore bi- 
ased) VLMO sample provides an upper limit for 
the 24 /im median flux of the VLMOs. If the ob- 
servations were not missing the fainter sources, the 
median 24 /im flux density would be lower. 

In order to quantitatively interpret the dif- 
ferences observed in the IRAC colors and in 
the distribution of the SEDs we used a grid 
of disk models with and without dust settling. 
We applied the RADMC radi ative transfer code 
( Dullemond fc Dominikll2004a ). which can deter- 
mine the vertical density structure of the disk 
self-consistently, assuming perfect dust-gas cou- 
pling and hydrostatic equilibrium for the gas. The 
RADMC is complemented with RAYTRACE, 
a post-processing tool to calculate the emerg- 
ing SED. The primary input parameters of the 
RADMC/RAYTRACE codes are the tempera- 
ture (T ff), luminosity (Lboi), the mass (M sta r) of 
the central star, and the inner radius (-Rin), outer 
radius (i?disk), inclination (i), and mass (Mdisk) of 
the disk. Another key parameter of the RADMC 
code is the number of photons used in the radia- 
tive transfer. The low photon number could cause 
fluctuations of the mid-plane temperature of the 
disk, especially at the inner 0.5 AU where the near 
and mid-infrared radiation of the dust is coming 
from. The fluctuation co uld be reduced b y a pho- 
ton diffusion algorithm ( Min et al.l 120091 ) , but it 



can not be used in every case. Thus, the photon 
number should be set carefully. The changing of 
the pressure scale height (H p ) with the radius de- 
scribes the disk flaring. The pressure scale height 
in RADMC is expressed in the unit of the distance 
from the star (H/a, where a is the radial distance 
from the star) and can be set externally. A flatter 
disk has a smaller H/a than a flared disk at the 
same emitting region. 

We adopted the median effective temperature 
and bolometric luminosit y of the s t ars of the tw o 
groups, as determined in iLuhmanl ( 20071 . 1200 8b). 
The stellar masses were derived from the theoret- 



the low-mass sample are shown in Table [2] 

First, we calculated flared disk models by set- 
ting the vertical structure self-consistently, i.e. as- 
suming perfect dust-gas coupling and hydrostatic 
equilibrium for the gas. The model grid param- 
eters and their range are as follows: R m from 2 

i?star tO 7 i? st ar, M disk from 5xl0~ 3 M sta r to 0.05 

M s tar- We kept the i?disk constant: in the case 
of the VLMO sample it was 50 AU, while in the 
case of the low-mass sample it was 100 AU. The 
pressure scale height was calculated by the code 
in 10 iteration steps. The number of photons was 
set to 300,000; this photon budget was found to be 
sufficient to achieve a smooth mid-plane temper- 
ature profile even in the inner regions of the disks, 
when using the diffusion algorithm. We calculated 
SEDs for inclinations ranging from 40° to 70° with 
5° steps. 

We performed a reduced \ 2 fitting for the me- 
dian SEDs to identify the best-fitting model. 
The errors were calculated by i^err.med. (A) = 
median(F c „(X)/F(X)) x median(F '(A)) . The 
best-fitting model for the low-mass group is 
R m = 2 J R star , M disk = 5 x 10~ 3 M st ar, and i = 55° 
with x 2 = 270.32. In the case of the VLMOs the 

best-fit is i? in = 2 J R star , Mdisk = 5 X 10" 3 M sta r, 

and i = 60° \ 2 = 205.57. The H/a pressure 
scale height of a fully flared, radiative equilibrium 
disk is given by H/a = v /fcTa 3 //ig as m p GM s t ar /a, 
where T is the mid-plane temperature at a radius 
of the disk, /x gas is the mean molecular weight 
of gas and m n is the photon mass (see eq. 6 in 
Dullemond fc Dominiktl2004bh . Thus, the H/a at 



ical evolutionary models of iBaraffe et al.l (J1998) 



a = 1 AU for the best-fitting models are 0.041 
(low-mass) and 0.085 (VLMOs). The best-fitting 
fully flared models are shown on Fig. [3] (upper and 
lower left panels). 

The fully flared models could not fit the median 
SEDs of the two groups satisfactorily within the 
parameter ranges of our model grid. The best- 
fitting model of the low-mass group overpredicts 
the median observed 24 /im fluxes, suggesting that 
a flatter disk is required. On the other hand the 
model underpredicts the SED at H, K s and IRAC 
wavelengths, this effect likely resulting from the 
absence of accretion in our models. To verify this 
we compared the amount of H and K s band flux 
excess to th e H an d K s fluxes of the accretion disk 



The adopted stellar parameters of the VLMOs and 



models of D'Ales sio et al.l ( 2004T) . The amount of 



excess at H and K s bands could be reproduced by 



Table 2 
Adopted stellar parameters. 


Parameters 


VLMO 


Low-mass 


r cff [K] 

ibol [Lq] 
M star [M®\ 


3,024 
0.03 
0.08 


3,669 
0.39 
0.60 



models with an accretion rate of ~ 10~ & M so \/yr, 
a typical mass accretion rate for a few Myr-old 
star with disk. Thus, the poor fit of the shorter 
infrared wavelengths is likely caused by the lack of 
accretion in our models. 

Fully flared disks provide an even worse fit for 
the VLMO group. Fitting the K s and IRAC data 
points requires unrealistically low disk mass and 
very extreme inclinations; on the other hand, these 
extreme models predict more flux at 24 /mi than 
the observed upper limit. 

In order to verify whether dust settling and flat- 
ter disks may explain the observed disk structures, 
we expanded our model grid with an additional pa- 
rameter. We express dust settling by allowing the 
pressure scale height at i?disk to vary between 0.01 
and 0.3. The flatter disks require more model pho- 
tons: we set the number of photons to 2.5 million 
to ensure continuous and smooth mid-plane tem- 
perature profile. On this expanded grid the best- 
fit model for the low-mass sample is i?; n = 2i? star , 
M disk = 0.02M star , % = 40°, and H/a = 0.05 with 
X 2 = 123.37. The best-fit model for the VLMOs 
is R in = 4i? star , M disk = 0.05M star , i = 40°, and 
H/a = 0.05 with x 2 = 5.89. Thus, allowing dust 
settling and flatter models provides much better 
fits and disk masses m ore consistent with those ob - 
served (see Fig. 8 in lAndrews fc Williams! 12005) 



and inclinations, close to what is expected from 
a large sample of randomly oriented disks. The 
best-fit models from this grid are shown in Fig. [3] 
(upper and lower right panels). 

Our expanded grid reveals that flatter disks are 
a much better match to the SEDs than flared ones. 
In the case of the low-mass group, the best-fit 
flatter disk model (Fig. [3] upper right panel) un- 
derestimates the H , K s , and IRAC fluxes similarly 
to the fully flared model, but it fits the 24 /im flux 
well. In the case of the VLMOs the best- fit model 



fits the observations almost perfectly. There is a 
little deviation at the IRAC [8] band. This non- 
perfect fit is expected, because the flux density at 
wavelengths close to 9 /zm are strongly affected by 
the shape of the silicate emission feature, a pa- 
rameter not varied in our models. The 24 /mi flux 
of the models is a sensitive proxy of the pressure 
scale height (and flaring) of the disk: a lower scale 
height results in a lower 24 /im flux. Thus, the 
upper limit of the 24 /im flux of the VLMOs trans- 
lates to an upper limit on the pressure scale height. 

To compare the pressure scale heights of the 
best-fit flatter models, we calculated the scale 
heights of the model disks at the radii where 
the observed IRAC and MIPS fluxes emerge. As 
a first approximation the protoplanetary disks 
could be described as a black body radiator with 
its temperature decreasing with the radius. Ac- 
cording to Wien's displacement law we calcu- 
lated the temperatures of the black bodies which 
have the maximum of radiation at the IRAC and 
MIPS wavelengths and interpolated the mid-plane 
temperature profile of the disk to determine the 
radius where the radiatio n emerges. Following 



Chiang fc Goldreichl (J1997I ) we assumed a flaring 



angle of 2/7 and calculated the pressure scale 
height at these radii (a) by H/a — [H/a]o x 
(^-) 2 / 7 , where ao and [H/a]o are the outer radius 
and the pressure scale height at the outer radius. 

Table [3] summarizes the pressure scale heights 
of the best-fitting fully flared low-mass, VLMO 
and flatter low-mass, VLMO disk models at the 
regions of the disks where the 3.6, 4.5, 5.8, 8, and 
24 /j,m radiation emerge. The scale height of the 
VLMO disks is reduced by about a factor of 4 
relative to their predicted structure. In contrast, 
that for the low-mass group is reduced by only a 
factor of 2. 

Our SED modeling also illustrates the power 
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Fig. 3. — Best-fit models with dust settling explain better the median SEDs than the best-fit fully flared 
models. The upper panels shows the best-fitting fully flared (upper left) and flatter (upper right) model 
of the low-mass group, the lower panels shows the best-fitting fully flared (lower left) and flatter (lower 
right) model of the VLMOs. Dashed lines represents the black body SEDs calculated with the adopted 
stellar parameters. The solid lines shows the best-fitting models. Note, that the relatively high \ 2 values 
for certain models are partly due to the non-black body spectrum of the stellar photosphere and the fine 
structure of the silicate emission feature, none of which is included in our models. 



Table 3 
Pressure scale heights of the best-fit disks at regions that dominate emission in various 

wavelengths 



Group 3.6 /mi 4.5 /im 



/jm 



8 /-im 



24 fim 







Fully flared models 






Low-mass 


0.0112 


0.0120 0.0128 


0.0135 


0.0178 


VLMO 


0.0181 


0.0188 0.0197 


0.0215 


0.0281 



Flatter models 

Low-mass 0.0046 0.0050 0.0054 0.0062 0.0088 

VLMO <0.0048 <0.0050 <0.0053 <0.0057 <0.0082 

Relative pressure scale heights 

Low-mass 2.42 2.43 2.37 2.18 2.02 

VLMO >3.79 >3.75 >3.73 >3.76 >3.42 



of mid- to far-infrared photometry in accurately 
determining the disk geometry and following the 
settling of dust and its impact on disk evolution. 
The recently launched Herschel Space Observatory 
provides a new opportunity to study disks at these 
hard-to-access wavelengths. In the following we 
point out that Herschel is capable of mapping the 
SED of typical VLMO disks in just minutes. We 
used our model SEDs to estimate the flux density 
at the 60-85 /im band of the Phot odetector Array 
Cam era and Spectrometer (PACS. IPoglitsch et al.l 
20081 ) of the Herschel Space Observatory. The pre- 
diction of flux density at this band for the best- 
fitting fully flared low-mass model is 1278 mJy, 
while for the fully flared VLMO model it is 117 
mJy. The best-fitting flatter low-mass model pre- 
dicts 424 mJy, while the best-fitting VLMO flatter 
model predicts 36 mJy. The instrument is capa- 
ble of 8 a detection in the faintest case with ~11 
minute exposure time. Thus, PACS measurements 
could distinguish between the fully flared and flat- 
ter geometries. 

5. Discussion 

Using colors derived from Spitzer photometry, 
median SEDs, and best-fit disk models we showed 
that the flaring of disks around young low-mass 
and very-low-mass stars (VLMOs) is statistically 
different and that disks around cooler stars are 
flatter. While previous studies by our group (e.g. 
Apai et al.ll2005t IPascucci et al.l 120091 ) and others 



already hinted at such a difference, its statisti- 
cal significance has not yet been verified. Dust 
settling - a physical process that leads to flatter 
disks - enhances solid surface density in the disk 
mid-plane and is an early, important step toward 
planet formation. 

Our study shows that in coeval disks this pro- 
cess is more efficient around lower-mass stars, an 
effect that has neither been predicted nor easily 
understood — in fact, theoretical works predict 
the opposite due t o the lower gravity of the cen- 



tral star (see, e.g. IWalker et al.ll2004l ). 



Although developing a detailed, quantitative 
model explaining the stellar-mass-dependent dust 
settling is beyond the scope of the current paper, 
we venture to speculate about possible explana- 
tions of the observations. Because the scale height 
of the fine dust is set by the dynamical coupling of 
dust and gas processes that lead to differences in 
the evolution of the gas disk, the dust disk or the 
dust-to-gas coupling may all lead to the observed 
stellar type-dependent disk structures. 

A lower gas surface density in disks around 
cooler stars may provide a straightforward expla- 
nation for lower dust equilibrium scale heights, i.e. 
flatter disks. This model would require that ei- 
ther disks around cool stars lose their gas faster or 
that lower-mass stars form from cloud cores with 
a lower-than-average gas to dust ratio. 

Other possible explanations include a weaker 



turbulence in disks around lower-mass stars, lead- 
ing to weaker coupling between the dust and 
gas components. As turbulence and viscosity 
also directly determines the average accretion rate 
through a disk, future studies may be able to find 
a consistent framework for connecting variations 
in disk structure and accretion with spectral type 
to systematic differences in disk viscosity. Because 
the origin of turbulence (and viscosity) in proto- 
planetary disks is poorly understood, observations 
such as those presented in this paper may provide 
indirect clues constraining on these long-standing 
questions. 

However, more complex explanations are also 
possible. Detailed simulations of grain-grain colli- 
sions and their impact on the overall disk structure 
and spectral energy distributions have been dis- 
cussed in the literature extensively, including the 
effec ts of collisions, fragmentation, inward drift 



(e.g. I Wcidcnschilling 1 997c iDominik fe Dullemond 
2008; Bi rnstiel et al1l2009h . Supp orted by labora 



tory studies of grain collision (e.g. I Blum fc Wurm 



2008) the models find a very efficient and rapid 



grain growth in disks on time scales shorter than 
the disk lifetime. Indeed, most models would 
lead to a rapid loss of the fine due component, 
in stark contrast to the observed large frequency 
of dust disks at even 1-3 Myr ages. This sug- 
gests that processes other than simple grain-grain 
collisions and fragmentation are involved in defin- 
ing the grain size distribution in disks (see also 
Pontoppidan fc BrearlevlbOlOl ). 



We also point out that the flaring disk surface 
is defined by a small mass fraction of the dust in 
the disk; and disks around the coolest stars and 
brown dwarfs will be about an order of magnitude 
less massive than those around sunlike stars. How- 
ever, the optical depths probed by infrared obser- 
vations will be independent of the stellar spectral 
type. In addition, if the dust component is in equi- 
librium then the scale height of the smaller dust 
component will be larger than that of the larger 
grains. Thus, the observations at the same wave- 
lengths may probe deeper in the less massive disks 
around cool stars and lead to a larger deduced typ- 
ical grain size and a flatter disk, even if the grain 
size distribution are identical in the two disks. 

The observed differences in disk structure 
is one of the first steps toward distinguishing 
between the interesting possibilities identified 



above. Even more importantly, however, these 
results question the assertion that the initial 
conditions for planet formation are not stel- 
lar mass-dependent, frequently made in current 
planet formation models. The observed differ- 
ence in the disk structure is a very important 
clue that joins th e ranks of differences observed 
in disk mass (e.g. Klein ct al. 20031; IScholz et al 



2006c IWilliams fc Jona than 2007), d ust pr oper- 
ties (lApai et all l2005t IPascucci et all 120091). ac- 



cretion rates (INatta et al.ll2004t iMuzerolle et al 
12005c iMohantvet al.l l2005h . Most recently, sta- 
tistically significant stellar-mass-dependent dif- 
ferences have been identified even in th e organic 
chemistry of disks (|Pascucci et al.ll2009l ). 

The evidence presented here on stellar-mass- 
dependent disk structures completes the picture 
and further emphasizes the need for disk mod- 
els that naturally reproduce these stellar-mass- 
dependent disk properties and disk evolution. 
Only with such models we will be able to develop 
predictive, quantitative models for planet forma- 
tion and understand the diversity of planetary 
systems. 

6. Summary 

Using Spitzer/IRAC colors of a complete sam- 
ple of coeval young stars in the Cha I star-forming 
regions we found that the IRAC colors of disks 
around the lower-mass and higher-mass stars are 
statistically different: disks around cooler stars 
are bluer, i.e. their SEDs have a less steep slope 
beyond 4 /im. The same difference is seen in 
the median SEDs of the two samples. Using a 
grid of radiative transfer disk models we found 
that: A) fully flared disks are not consistent 
with the median SEDs; B) disks in the lower- 
mass sample require a much stronger reduction of 
the disk scale height relative to the fully flared 
disk geometry than those in the higher-mass sam- 
ple. In short, the Spitzer IRAC data reveals a 
stellar-mass-dependent disk structure and iden- 
tifies more efficient dust settling in disks around 
lower-mass stars than in their higher-mass coun- 
terparts. These findings demonstrate that the ini- 
tial and boundary conditions of planet formation 
around cool stars are different from those around 
sun-like stars. 
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